arXiv: 1507.04514v2 [cond-mat.mes-hall] 6 Jan 2016 


Decoherence of a single spin coupled to an interacting spin bath 
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Decoherence of a central spin coupled to an interacting spin bath via inhomogeneous Heisen¬ 
berg coupling is studied by two different approaches, namely an exact equations of motion (EOMs) 
method and a Chebyshev expansion technique (GET). By assuming a wheel topology of the bath 
spins with uniform nearest-neighbor AA-type intrabath coupling, we examine the central spin dy¬ 
namics with the bath prepared in two different types of bath initial conditions. For fully polarized 
baths in strong magnetic fields, the polarization dynamics of the central spin exhibits a collapse- 
revival behavior in the intermediate-time regime. Under an antiferromagnetic bath initial condition, 
the two methods give excellently consistent central spin decoherence dynamics for finite-size baths 
of A < 14 bath spins. The decoherence factor is found to drop off abruptly on a short time scale and 
approach a finite plateau value which depends on the intrabath coupling strength non-monotonically. 
In the ultrastrong intrabath coupling regime, the plateau values show an oscillatory behavior de¬ 
pending on whether A/2 is even or odd. The observed results are interpreted qualitatively within 
the framework of the EOM and perturbation analysis. The effects of anisotropic spin-bath coupling 
and inhomogeneous intrabath bath couplings are briefly discussed. Possible experimental realization 
of the model in a modified quantum corral setup is suggested. 

PACS numbers: 03.65.Yz, 72.25.Rb, 73.21.La 


I. INTRODUCTION 

A single central spin coupled to a spin environment [1] 
is a widely studied quantum many-body system because 
of its high relevance to solid-state based quantum com¬ 
putation [2-32], quantum decoherence and quantum in¬ 
formation [33-42], and excitation energy transfer [43-45], 
etc. In these systems, a single spin is inevitably coupled 
to its surrounding (pseudo-)spin environment and it is 
important to understand and control the decoherence of 
the central spin induced by environmental fluctuations. 
Among these, the central spin or Gaudin model [46] plays 
an important role in describing the decoherence dynamics 
of an electron trapped in a quantum dot, which has long 
been proposed to be a promising candidate for realizing a 
qubit for solid state quantum computation [47, 48]. The 
trapped central spin S coupled to a bath of N nuclear 
spins {Sj} via inhomogeneous hyperfine couplings {gj} 
is well modeled by the Hamiltonian 

N 

H = (jjSz + ■ Sj, (1) 

i=i 

where w = gg-sB is the Larmor frequency of the central 
spin due to the externally applied magnetic field B. 

The dynamics of the central spin model Eq. (1) has 
been studied theoretically with a wide variety of meth¬ 
ods. In the strong magnetic field limit u) > NiniBx{gj}, 
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perturbative treatments [2-4], non-Markovian master 
equation methods [5-8], and effective Hamiltonian ap¬ 
proaches [9-11] in terms of the flip-flop hyperfine interac¬ 
tion 9j{S+S~ + S-S ^) were commonly employed. In 
the low-field limit, semiclassical approaches based on the 
spin coherent state representation [12-14], exact meth¬ 
ods via either the Bethe ansatz techniques [15-18] or 
exact diagonalization [19-21] have been used intensively 
to study the dynamics of the model. Most recent ad¬ 
vances include a combination of algebraic Bethe ansatz 
and Monte Carlo simulation [22, 23], as well as an ad¬ 
vanced numerical method based on the Chebyshev poly¬ 
nomial technique [24]. The extension of the density ma¬ 
trix renormalization group [49] to the real-time dynamics 
in central spin model [25] can treat up to A « 1000 nu¬ 
clear spin, but is restricted to rather short times. 

In the above studies, the dipolar interaction between 
nuclear spins was often neglected because it starts to play 
a role only on longer time scales compared to the hyper¬ 
fine interaction. However, these effects need to be taken 
into account if one is interested in the long-time behav¬ 
ior of the central spin dynamics. On the other hand, it 
is also of interest to study the influence of intrabath in¬ 
teractions [31-36] on the decoherence properties of the 
central spin from a purely theoretical point of view. The 
bath spins described by Eq. (1) form a noninteracting 
spin cluster whose spatial correlations are only reflected 
from its coupling {gj} to the central spin. 

In this work, we use exact equations of motion ap¬ 
proach combined with a Chebyshev expansion technique 
to study the decoherence dynamics of a central spin cou¬ 
pled to a periodic A A spin chain with nearest-neighbor 
interaction. Such a system can be obtained by arranging 
the bath spins uniformly on a ring but leaving the hyper- 
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fine interactions unchanged, then by turning on the ho¬ 
mogeneous XX-type interaction between adjacent bath 
spins. We note that such a fictitious one dimensional 
structure of the spin bath also helped implement density 
matrix renormalization group algorithms of Eq. (1) [25]. 
The present system can also be thought of as an extension 
to the model studied in our previous work [36], where we 
obtained the exact reduced dynamics of a central spin 
coupled to an XX spin ring via the XX-type interac¬ 
tions as well. Taking advantage of the translational in¬ 
variance of the spin bath as well as the conservation of 
the total magnetization, we can write down the evolved 
state of the system in terms of Jordan-Wigner fermions 
in the momentum space of the chain. Interesting quan¬ 
tities like the decoherence factor and polarization of the 
central spin then can be expressed as combinations of the 
corresponding time-dependent probability amplitudes. 

We have studied the spin decay of an initially fully po¬ 
larized central spin as well as the decoherence factor [50]. 
For the central spin model described by Eq. (1), due to 
the integrability and the conservation of total magneti¬ 
zation in z-direction, a partial coherence is maintained 
[22-25, 31] depending on the initial preparations of the 
system [22, 23]. Switching on the XX-ring coupling de¬ 
creases the long-time limit of the decoherence factor and, 
therefore, the partial coherence in the system. Once the 
XX-ring coupling exceeds some value, however, we find a 
gradual increase of coherence, which appears to be coun¬ 
terintuitive at first. We present a physical picture for 
explaining these surprising findings. 

Although the EOM method employed here restrict us 
to small baths with N < 14, it provides benchmarks for 
testing the GET method, which can be applied to large 
baths. It has been demonstrated [24] that the GET ap¬ 
proach [51-55] is well suited to access the quantitive fea¬ 
ture of the central spin dynamics. In particular, quantum 
corrections [5, 26] to a quasi-classical treatment of the nu¬ 
clear spin bath [28] have been revealed. Furthermore, the 
EOM approach also provides a possibly intuitive way to 
understand the observed nontrivial dependence of the de¬ 
coherence on the intrabath coupling strengths. We finally 
discuss potential reality of the design in the laboratory. 

The rest of the paper is structured as follows. In Sec. 
II we introduce our model Hamiltonian and describe the 
EOM approach and the GET technique in detail. In 
Sec. Ill we study the central spin polarization for a 
fully polarized bath, and the decoherence dynamics for a 
zero-polarization bath initial state by employing the two 
methods. Gonclusions are drawn in Sec. IV. 


II. METHODOLOGY 


type interaction between adjacent bath spins. Such a 
system can be described by the Hamiltonian 

H = Hs + Hb + Hsb, 

Hs = ujSz, 

N 

HB = jJ2{s]sj+, + spy^,), 

N 

Hsb = + SyS]) + g'S^S^l ( 2 ) 

i=i 

where we have separated H into the system part (cen¬ 
tral spin) Hs, the environment part (XX ring) Hb, and 
the anisotropic Heisenberg system-bath coupling Hsb be¬ 
tween the two. In Eq. (2), Sa and Sf are the spin-1/2 
operators for the central spin and spin j in the XX bath, 
respectively, uj is the Larmor frequency in an external 
magnetic field, J is the uniform nearest-neighbor intra¬ 
bath coupling strength, and {gj} as well as {g'} are the 
exchange interactions between the central spin and the 
bath spins. 

The Hamiltonian (2) includes the anisotropic spin 
model [26, 27] which corresponds to A = g'j/gj and is rel¬ 
evant to quantum dots charged with a single hole, where 
A depends on the mixing of heavy and light holes. For 
A = 1, the isotropic model describing electron charged 
quantum dots is recovered. The Hamiltonian (2) also ex¬ 
tends the model proposed in Ref. [36] by including the 
Ising component of the system-bath coupling, which is 
essential for the model to cover the Gaudin model in the 
limit of J —^ 0. It is thus experimentally more relevant 
to introduce the Ising part of the system-bath coupling. 
As we will see, the decay of the central spin decoherence 
can be suppressed by increasing A. 

It can be easily seen that the total magnetization 
M = Sz + Lz is conserved, where Lz = is 

magnetization of the spin bath. Thus, we dropped the 
Zeeman term for the bath spins in Eq. (2). However, 
the model described Eq. (2) is only integrable for J = 0 
and A = 1, i.e., for the conventional Gaudin model with 
a noninteracting spin bath and isotropic spin-bath cou¬ 
pling. Hence, results from the Bethe ansatz solution gen¬ 
erally cease to be applicable here. 

In order to separate dynamics governed by the fluc¬ 
tuations of the Overhauser field [28] from the dynamics 
induced by the spin-bath interaction, we define the en¬ 
ergy scale Wfluc [24] 


^fluc 


N 




( 3 ) 


A. The interacting central spin model 

By arranging the N spins (X even) in the spin bath 
into a ring-like fashion, we consider an extension of model 
(1) by introducing a homogeneous nearest-neighbor XX- 


associated by the fluctuation of the Overhauser field. The 
such measured short-time dynamics will be universal for 
different number of bath spins, which can be physically 
understood from second-order short-time perturbation of 
the dynamics. Regardless of the initial conditions, the 
second-order process flipping the same spin back always 
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yields a dimensionless term oc Below, we 

will measure all energies in units or Wfluc- 


B. The exact equations of motion approach 


The one-dimensional bath Hamiltonian can be di¬ 
agonalized by using the Jordan-Wigner transformation 
S'- = S? -iSf 


n^(l - 2 cJcj)c,;, 


Sf = C^C,- — i 


with c| a fermion creation operator. Then H describes a 
central spin immersed in a spinless fermion bath 


H = uj'S^ + + 4+i^o) 


N 


j=l 


N 


- fe ic]T,S. + c,t]S+) + 25 'S,c]c,], (4) 
j=i 


where the Jordan-Wigner string Tj = exp(i 7 r^;^j^ cjci) 
and 


N 

J =u:-Y,9'j 

j=i 


( 5 ) 


is the effective magnetic field experienced by the central 
spin when all the bath spins are pointing down. De¬ 
pending on the parity of the eigenvalue Nf oi the total 
fermion number operator A// two subpaces 

with projection operators P+ = even Nf) 

and P- = identified. This 

results in antiperiodic cjv+i = —ci or periodic boundary 
conditions cn+i = ci imposed on the fermions for even or 
odd Nf. As a result, two sets of Fourier transformations 
can be introduced 

^ E («) 

k&K± 


where {c±^k} are Fourier modes with wave numbers sur¬ 
viving in k± = {-TT-h (4j - 3)^ ± 2 ^}, j = l,...,iV. 
Now Hb is diagonalized as 


Hb=Y1 

fT—it 

Hu = ^ ^ (7) 

keK„ 

where Sk = J cos k is the single-particle spectrum. We 
will work in the interaction picture with respect to the 
first two terms of Eq. (4). Direct calculation gives 


N 

Hi{t) = Y. PuP^-k]Tje-^^--^P-u 

j=l u=± 

N 


-fh.c.) J- 2 


i=i 
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Pu. 
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We assume that initially the central spin and the bath 
are decoupled (the up and down states of each spin are 
denoted by | 1 ) and | 1 ), respectively) 

l'0(O)) = (ai|i) + ai|l))l<(>fe), (9) 

where the coefficients satisfy |aip -I- |aip = 1 and \(f)b) 
is the bath initial state which is assumed to have fixed 
magnetization = m — ^ with m the number of up 
spins in the bath. For \(j)b) having nonzero overlap in 
several L^-subspaces the time evolution obtained in dis¬ 
tinct subspaces can be superimposed. 

It is convenient to define two vectors made up of mo¬ 
mentum and spatial indices 

km = (^11 ^2 ; ■•■5 km )^ Jm = (jl; J2 5 ^771)5 ( 19 ) 

with the convention ki < k 2 < ■■■ < km and 1 < ji < 
j 2 < ■■■ < jm < N, and the total energy of the m fermions 
occupying the set of modes km 

m 

^km ~ E 
1=1 

Then \(j)b) can be written in the momentum space 

\4>b) = ^k^ I^c7,7n); ( 12 ) 

km 

where \k^^m) = |9) with a = -l-(—) for even 

(odd) m. Here |0) is the vacuum state of the fermions 
corresponding to the fully polarized bath state | 1 ...I). 
Normalization of jt/fb) gives 

EI^S7„I' = 1- (13) 

km 

Since the total magnetization is conserved, the time 
evolved state in the interaction picture is of the form 

\^iit)) = 

km 

fem + l 

+ E aiC,-^_,(l)|l)e”^^--^‘|fc-.,7n-i). (14) 

fem-l 

The nonvanishing initial values of A, B, C, and D can 
be read from Eq. (12) as 

^km^^)=Bkm^^)=^km 

Applying the Schrddinger equation f9t|'0/(t)) = 

Hi{t)\ipj{t)) and after some algebra we arrive at the fol¬ 
lowing sets of EOMs for the probability amplitudes A, 
B, C, and D 

iC'pm-l = ^Pm-lP'Pm-l + (E^fem® fkm-,Pm-l 

km 

km — 1 


( 8 ) 


(16) 
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'^^Pm — ^Pm^Pm + ( ^km-1^ ^Pm;km-1 

km-1 


and 


*-®Pm — %m^Pm + ( ^ ^k^ + 1^ ^k^ + lWr, 


fcm + 1 


'^^ k ^ hm -, Prr )^ 


(18) 


iDpm + 1 — ^Pm + l^Pm + 1 + ^fc, 

km 

- £>- /■- ) 
fcm+l''fcm+i;pm + i'' 


g-iw't 

Pm + l j^TT 


(19) 


The derivation of these EOMs for the case of g' = 0 can 
be found in Ref. [36] and extension to the present model 
is straightforward. The auxiliary functions / and f here 
are defined to be 

n+1 

fkn + i-,pn(^9j}) — ^ ^ ^kn + i;jn + i 

Jn + 1 






931 ^.^(0 ! ( 20 ) 

PnJri + 1 


and 


( 21 ) 


where the Slater determinant .j = det(O) with the 
n X n matrix O made up of plane waves 


Oa,b = (a,&= l,2,...,n) 

\/N 


( 22 ) 


The vector j^^ri = ih, in Eq. (20) 

is a string of length n with the element ji being removed 
from j„+i. 

We now turn to some remarks regarding the above pro¬ 
cedures. (i) We note that our treatment here is similar to 
that in Ref. [30] where the real-space basis states of the 
bath were used to study the dynamics of model (1); (ii) 
The system-bath coupling configurations {gj} {{g'j}) are 
incorporated into the /-functions (/'-functions), while 
the intrabath coupling J enters the equations of motion 
through the occupation energy . Eurthermore, the 
time-dependent phase factors * associated with the 
/-functions depend both on the bare magnetic field uj and 
average of the transverse couplings {g'j}] (iii) The num¬ 
ber of / or /' functions is of the order of ~ (C'jv)^ (here 
= {N-my. binomial coefficient), as a result, it 

is usually time consuming to numerically calculate these 
functions for m ~ N/2] (iv) Intriguingly, the same struc¬ 
ture with the /-functions in Eq. (20) also appeared in the 


study of nonlinear optical response of one-dimensional 
molecular aggregates [56]. 

We are interested in the polarization {Sz{t)) and the 
decoherence factor \r{t)\ = |(S'+(t))/(5'+(0))| [50], which 
can be easily calculated by transforming Eq. (14) back 
to the Schrbdinger picture 

(s.w)= 

km—l km 

+ 5 l“-t(El%J®- E |OE.„t). ( 23 ) 

km km+1 

h^)I = iE4^.(^)^LWI' (24) 

km 


Using Eq. (14), it can be checked that {Sz{t)) is two times 
the autocorrelation function [3, 25] 

^(t) = (^(0)|5,(t)5,|^(0)) = i(5,(t)). (25) 

Furthermore, as a direct consequence of the EOMs 
Eqs. (16)-(19), one can derive the following relation 


d\r{t)\ 

dt 


|t=o 


= 0 , 


(26) 


regardless of the bath initial conditions. Eq. (26) states 
that the short-time decay of |r(t)| begins at least with 
the quadratic term of t, e.g., a Gaussian initial decay of 
|r(t)| was observed in Ref. [36]. 

Before ending this section, let us see how the EOMs 
look like for a uniform hyperfine coupling. In this special 
case, we have gj = g and g} = g', then the /- and /'- 
functions Eqs. (20) and (21) can be evaluated in closed 
form [57] 


n+1 n 

(9) = (27) 

j=l i=l 

,, ) n,->/ ) 

nLin;;i(i-e-®('= 2 -pd) 

and 

n 

(28) 

i=l 

Eqs. (27) and (28) show that the spin-flip part of the 
anisotropic hyperfine coupling will drive the initial state 
l+(0)) = Efe„(ai|i)+ai|l))Xfc„l+m) into states within 
the (m± l)-sectors that conserve the total initial momen¬ 
tum k, while the Ising component only induces 

transitions into the states within the same m-sector that 
preserve the exact momentum configuration km- This 
arises from the circular symmetry of the interaction con¬ 
figuration for the uniform hyperfine coupling. For general 
nonuniform coupling, however, no analytical results for 
the /- and /'-functions exist and we have to resort to 
numerical methods. 
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C. The Chebyshev expansion technique 

The Chebyshev expansion technique (GET) [51-54] 
has been developed for more than three decades and of¬ 
fers an accurate way to calculate the time evolution of a 
general single initial state |' 0 o) under the influence of a 
general time-independent Hamiltonian H with a bound 
spectrum: 

|^(<))=e-*^‘|^o). (29) 


Hamiltonian H' . Since Jn{x) ^ {exj2n)'^ for large or¬ 
der n, the Chebyshev expansion converges quickly as n 
exceeds AEt. This allows to terminate the series (33) 
after a finite number of elements Nq guaranteeing an ex¬ 
act result up to a well defined order. More details on 
the method can be found in Refs. [24, 53, 55] and in the 
review on kernel polynomial methods [54]. 

III. RESULTS 


The main idea of the method is to construct a stable 
numerical approximation for the time-evolution operator 
jg independent of the initial state \tpo), whose 
error can be reduced to machine precision for any given 
time t. While this can be achieved by any complete set of 
orthogonal functions, the Chebyshev polynomials guar¬ 
antee a linear dependence of the time t and the order of 
the expansion. 

If the spectrum of the Hamiltonian is bound to Emin < 
E < Emax, the time-evolution operator can be 

expanded as 

OO 

^-^m ^ bn{t)T^{H') (30) 

n—0 


A. Fully polarized bath: application of the EOM 
approach 

As an illustration for the EOM approach, we first ap¬ 
ply it to the simplest initial bath configuration, the fully 
polarized bath [2, 3, 16], 

N 

]</),) = 11...1) = n4l0) = (35) 

1=1 

The central spin is assumed to be in its down state jl) 
initially, i.e., oj = 1 and ai = 0, so that only A and 
the N C’s are used. Eq. (35) yields the following initial 
conditions for the coefficients 


after mapping the Hamiltonian to the dimensionless 
H' = {H — a)/AE where we have defined the center 
of the energy spectrum a = (Emax + Emin)/2 and its 
half-width AE = (Emax — Emin)/2. The time depen¬ 
dency is entering only the analytically known expansion 
coefficients 

b„(t) = (2 - J„(AEt), (31) 

where Jn{x) is the Bessel function of n-th order, while 
the dynamics of the systems is included in the Chebyshev 
polynomials r„(z) defined by the recursion relation 

r„+i(z) = 2zTn{z) - Tn-l{z), (32) 

subject to the initial conditions Tq[z) = 1 and Ti{z) = 
z. By applying the expansion of time-evolution operator 
(30) to the initial state j'0o) one obtains 

OO 

li’it)) = '^bn{t)\cl)n), (33) 

n—0 

where the infinite set of states 1^„) = T„(E')j'!/:o) obey 
the recursion relation[54] 

l<^„+l) = 2E'l(/)n) - |</>n-l), (34) 

subject to the initial condition j())o) = {ipo) and \4>i) = 

E'lV’o). 

The Chebyshev recursion relation of Eq. (34) reveals 
the iterative nature of the calculations. Starting from 
the initial state one constructs all subsequent states 
\4>n) using repeated applications of the dimensionless 


=SZ . , , = 0. 

PN'Jn’ PN-1 


(36) 


For the fully polarized bath, the dimension of the Hilbert 
space involved is A -|- 1 which allows for relatively large 
number of bath spins. To compare our results with prior 
works, we follow Ref. [16] to choose the following distri¬ 
bution of {gj} 

9j = j = 1,...,N (37) 


where g defines an overall energy scale. 

Let us first focus on the noninteracting and isotropic 
case with J/wfluc = 0 and A = 1. Fig. 1(a) shows the 
long-time dynamics of the polarization {Sz{t)) in the 
absence of the magnetic field, w/wfluc = 0 . 0 , obtained 
by solving Eq. (16) and Eq. (17) with initial conditions 
Eq. (36). The same result can be simply obtained by 
numerically integrating the equations of motion of real- 
space wavefunctions [2, 3]. The inset of Fig. 1(a) shows 
the initial evolution, reproducing the result in Ref. [16] 
through numerically solving the Bethe ansatz equations 
for model (1). For this case, we see that {Sz{t)) oscil¬ 
lates with a decreasing amplitude as time evolves, and 
the overall envelope finally decays in the long time limit. 
This is qualitatively consistent with Ref. [2, 3]. Our re¬ 
sults indicates that the analytic approximations derived 
from Bethe ansatz solutions [16] can correctly capture 
the polarization dynamics in the short or intermediate 
time regime, but may break down at sufficiently long 
time scales. 

For negative magnetic fields, a resonant field around 
which exhibiting maximal oscillation amplitude and 
minimal frequency can be defined through tOr = 


6 




FIG. 1. Polarization dynamics {Sz{t)) in the absence of intra¬ 
bath interactions, J/tonuc = 0. The bath is initially prepared 
in the fully polarized state Eq. (35) and the hyperfine inter¬ 
action is chosen to be isotropic with A = 1. The number 
of bath spins is = 30. (a) In the absence of the mag¬ 

netic field, cu/oifluc = 0. The inset shows the initial evolution 
up to Wfluci = 120. (b) LJ = 0.5 oJr, oJr, and 1.2 utr, where 
utr = —2.1353 cufluc is the resonant field dehned in Ref. [16]. 
(c) For positive magnetic fields oj/utfiuc = 0.5, 1.5, and 2.5, 
the polarization dynamics shows collapse-revival behaviors. 


(- Y.j ffi)wfluc [16]. In Fig. 1(b) we plot (S'z(t)) for three 
different negative magnetic fields uj = 0.5 uJr, and 

I. 2 UJr- We see that consistent results with Ref. [16] 
are obtained. Interesting dynamics happens for positive 
magnetic fields. Fig. 1(c) shows the polarization dynam¬ 
ics for w/wfluc = 0.5, 1.5, and 2.5. As the magnetic field 
increases, the Zeeman term loSz dominates and both the 
amplitudes of initial oscillations and the derivations of 
the overall profile from {Sz{0)) = —1/2 decrease. In¬ 
terestingly, collapse and revival behaviors [58] arise for 
intermediate and strong magnetic fields that are larger 
than w/wfluc ~ 0.6, during essentially the same period of 
time Wfluci G [0,700]. No obvious collapse and revivals 
occur after ajuuct ~ 700 and the dynamics become less 
regular. 

We also study effects of finite intrabath coupling J and 
different anisotropic factors A on the polarization dynam¬ 
ics. Fig. 2(a) presents {Sz{t)) for vanishing A, for which 
the model described by Eq. (2) is reduced to the model 
studied in Ref. [36]. Only results for positive J are shown, 
since the dynamics is symmetric under the sign change of 

J. Fig. 2(b) and (c) show the polarization dynamics for 
A = 0.5 and A = 1.0, respectively. To compare results 
for different A’s, we use the same energy scale Wfluc for 
A 7 ^ 1. We observe that: i) for fixed J, the oscillation 


FIG. 2. Polarization dynamics (S'z(t)) for different intrabath 
coupling J and anisotropic factor A, with the bath initially 
prepared in the fully polarized state Eq. (35) and in the ab¬ 
sence of the magnetic field, Lj/uJauc = 0. The number of bath 
spins is A = 30. (a) A = 0.0 with J/ujuuc = 0.0, 0.5, 1.0, and 
3.0, (b) A = 0.5 with J/ojauc = 0.0, ±1.0, and ±3.0, (c) A = 1 
with J/oJuuc = 0.0, ±1.0 and ±3.0. 


amplitude decreases, while the frequency increases as we 
increase the A. This means that the Ising component of 
the hyperfine interaction tends to increase the effective 
magnetic field on the central spin generated by the bath 
polarization, ii) for fixed A, amplitude increasing and fre¬ 
quency decreasing also occur for increasing J to positive 
values. However, for large enough J and finite A, fast os¬ 
cillating of {Sz{t)) tends to establish an envelope slowly 
oscillating around a mean value, iii) In the case of A 0, 
changing the sign of J to negative values will increase the 
overall amplitude of the oscillation profile to regions of 
{Sz{t)) > 0. 


B. Zero-polarized antiferromagnetic bath state 


We now turn to the study of the central spin deco¬ 
herence dynamics with the bath initially prepared in a 
zero-polarization antiferromagnetic state [22], 

l« = |iii-iii) = Es.T<2,4.(38) 

fejv ^ 

2 

yielding the initial values for the y's 


Xk, 


_ 

~ ^kN ;(2,4,...,Ar)‘ 

T 


(39) 


For this type of initial state, we will use an exponential 
distribution of coupling constants [22] 

q _ 


(40) 
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FIG. 3. (a) Real-time evolution of the decoherence factor 

|r(t)| starting from the zero-polarization bath initial state 
Eq. (38) for cu = 0, = 12 bath spins using the EOM 

approach. The intrabath coupling strength J/wauc are cho¬ 
sen to cover both the weak coupling and strong coupling 
regime. The dashed and dotted curves present the corre¬ 
sponding decoherence dynamics for A = 0.5 and A = 0.0, 
respectively. The inset shows the magnification of the in¬ 
terval UflucTmax € [100,200]. (b) |r(t)| for the same initial 

conditions but obtained using the GET approach and = 16 
bath spins, in the isotropic case A = 1. (c) The long-time 
(cjflucTmax = 700) mean value R{J) defined in Eq. (42) versus 
T/tufiuc for N = 8 and 16, in the isotropic case A = 1. 


where g defines the overall energy scale. In this case, 
Wfluc is given by 


^fluc — jY ^ 


„4r-l 



(41) 


The nonperturbative crossover regime and strong mag¬ 
netic field regime are defined by Wfluc ^ ^ ^ 2^ and 
w > 2g, respectively [22]. 

Here the main quantity of interest is the decoherence 
factor |r(t)| of a system initially prepared in a superposi¬ 
tion initial state ai|l) + aijl) for the central spin. Con¬ 
servation of the total magnetization restricts the bath 
polarization to the three sectors with m = y ~ 
and Y + 1; independent of J. The large number of /- 

functions (~ (C^'^^)^) to be calculated prevents us from 
dealing with baths even with intermediate number of 
spins. Therefore, we limit ourselves to A < 14 in the fol¬ 
lowing numerical calculations using the EOM approach. 

We have performed benchmark calculations in small 
systems by simulating the decoherence dynamics of the 
central spin using the two different methods. Fig. 3(a) 
and (b) show the time evolution of the decoherence factor 
|r(t)| for several different values of the intrabath coupling 
J at zero magnetic field, w = 0, with the number of bath 
spins A = 12 and 16, respectively. Since the EOM and 
GET results are identical for A = 12, Fig. 3(a) presents 
the data for both methods while we have employed only 
the GET approach for the results in Fig. 3(b) for A = 16. 
The intrabath coupling strength J are chosen to cover 
both the weak and strong coupling regime, by noting 
that the largest hyperfine coupling (in the isotropic case) 
2g/N Ri 0.42 Wfluc and 0.37 Wfluc for A = 12 and 16, 
respectively. 

In Fig. 3(a), we also present results for anisotropic 
hyperfine couplings with A = 0.5 and A = 0.0. We see 
that, for relatively small J, increasing of the Ising compo¬ 
nent of the hyperfine coupling strength tends to enhance 
the coherence. We also calculated |r(t)| for negative J’s 
(not shown), and find no essential deviation of the overall 
profiles from that of positive J’s. This might be due to 
the fact that an overall sign change of the bath Hamilto¬ 
nian does not change the energy gaps of the bath. In the 
following discussion, we will focus on the isotropic case, 
A = 1. 

The effect of a non-vanishing J 0 is significant. For 
all values of J considered, |r'(t)| drops off from the initial 
value of |r(0)| = 1 on the characteristic short time scale 
1 /wfluc set by the fluctuation of the Overhauser field and 
tends to oscillate around a fixed average value. Gonsis- 
tency with Ref. [22] is found for J = 0, where |r(t)| os¬ 
cillates about a mean value of « 0.55 at immediate time 
scales for A = 12. With increasing A, this mean value 
(defined below in Eq. (42)) gradually decreases as can be 
seen by comparing the J = 0 cases in Fig. 3(a) and 3(b), 
as well as from Fig. 3(c). We expect that it will reach 
1/2 as A —>■ oo if the long-time average of the coefficients 
|H| and |i?| remain finite and the initial spectral weight 
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Switching on the bath interaction J introduces addi¬ 
tional spin-flip processes in the spin bath. After such 
a spin-flip process occurs between two bath spins, how¬ 
ever, the Ising contribution of the hyperfine interaction 
will change sign leading to a destructive interference of 
the terms in Eq. (24). Consequently, the long-time limit 
of |r(f)| decreases with increasing J starting from J = 0 
as can be observed in Fig. 3(a) and (b). 

The other extreme regime of our model is the limit 
of very large J. In this case, the initial state of the 
bath \(j)b) is expanded in the eigenstates of the XX ring, 
\ka-,m) with fixed m = N/2, and the hyperfine coupling 
to the central spin is considered as a weak perturba¬ 
tion. For a finite the spectrum is discrete and the 
differences between different eigenenergies are increasing 
with J. Above some value of J, even the smallest ex¬ 
citation energy exceeds the largest energy scales of hy¬ 
perfine coupling, gi, and spin-flip processes correspond¬ 
ing to a creation of a particle or hole excitation will be 
suppressed. Consequently, the long-time limit of |r(t)| 
increases. For a system with no degeneracies, full coher¬ 
ence, i. e. |r(t)| = 1, would be restored, and quantum 
coherence is maintained at all times. However, spin-flip 
processes can couple energetically degenerate eigenstates 
with m and mil. One can Hnd different combinations of 
km and km±i such that their corresponding occupation 
energies are the same. However, since the dispersion for 
Nf even and odd are different, this coincident energy can 
only be = 0. Independent of the strength 

of the perturbation, gj j J, the degeneracy is lifted upon 
switching on the hyperfine coupling, and the states are 
maximally entangled. The arbitrarily weak spin-flip pro¬ 
cesses generated a mixture of local spin up and down 
states, and, therefore, the mean values of |r(t)| reaches a 
plateau less than unity for J —>■ oo. This plateau value 
is dependent on N. One may naively think that this 
plateau value will decrease with increasing N due to the 
increasing degree of degeneracies of coupled subsectors. 
However, our numerical simulation indicates that this is 
not the case, as we show below. 

To this end, the decoherence is increasing with increas¬ 
ing J/wfluc starting from J = 0, and the coherence is 
partially restored in a finite size system, for J —>■ oo due 
to the entanglement in energetically degenerate subsec¬ 
tors. Consequently, there must be a value J/uJduc where 
the long-time mean value of the decoherence factor |r(t)| 
approaches a minimum. In order to quantify this point 
we have defined the time averaged R{t) 

= - / dT\r{T,J)\ (42) 

t Emm Jtmin 

for t > tmin in order to eliminate the finite-size oscil¬ 
lations of |r(t)| in the long-time limit for different J. 
To separate the initial short-time decay of |r(f)| from 
its long-time mean values, we start the integration at 

^fluc^min — 100. 



FIG. 4. Plateau value of R{J) as function of N for J/ujauc 1 
calculated at a fixed J/uJuuc = 10 ®tt. 


We define the J dependent long-time limit as R{J) = 
I?(Tmax, J)- The results for R{J) at the largest time scale 
of our simulation, WflucTmax = 700, are depicted in Fig. 
3(c) as function of J/wfluc for A^ = 8 and 16. For J = 0 we 
note the predicted decrease of R{J = 0) with increasing 
N towards 1/2. Surprisingly, for N = 8, R{J) vs J shows 
not only a single, but several very sharp local minima: 
the first two at J/ujuuc = 0.28 and 0.47. This is a clear 
indicator of finite-size effects: at some values of J the 
decoherence is enhanced, which might be related to an 
increase of the degrees of degeneracy. 

For N = 16, R{J) does not exhibit such pronounced lo¬ 
cal minima, almost vanishes at J/ujuuc ~ 0.23 and slowly 
increases beyond that point. The number of eigenstates 
of the system are exponentially larger than for N = 8, 
and, therefore, the mean level spacing of the excitation 
spectrum is considerably smaller. Nevertheless, R{J) in¬ 
creases again for larger J after passing the minimum. 

We always find that R{J) approaches a finite mean 
value for large J. In Fig. 4, R( J) is plotted as function 
of A^ = 4, 6,8,10, and 12 at J/wfluc = 10®7r 1. We find 

an oscillatory behavior with N/2 changing its parity. In 
order to qualitatively understand this behavior, we note 
that there are totally = C// eigenstates in the m- 

sector of the bare XX ring. Among these, the number 
of states with zero energy is denoted by Then the 

spin-flip induced decoherence can be roughly measured 
by the following degeneracy factor 


0 


(™) = 
N — 


at;™) 




(43) 


The larger the large the decoeherence. Although 

is a regular function of N, the number of degen¬ 
erate states varies with N less regularly, espe¬ 

cially for odd N/2, due to different spectrum structures 
for N/2 = even or odd. For example, A4g^\ and 

AA^iJ all equal to 2, but Adjg^ = 692. Calculating 0^'^^^ 
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for N = 6 , 8 ,10, and 12 we get = 3.33 x 10 = 

4.59 X 10-2, ©Jq^ = 9.83 x 10-^, and ©J®^ = 7.62 x lO-^, 
which basically reveal the oscillating behavior of R{J) 
with increasing N. The details, however, will depend 
on the matrix elements. Since we recover i?( J) = 1 for 
gj = 0 , the decoherence is discontinuous for gj = 0 to 
finite gj > 0 which is related to lifting of degeneracies 
in the interacting spin bath and maximizing the entan¬ 
glement in this degenerate subsection by switching on gj 
independent of its values. 

Although the degeneracy argument presented here 
rests on a translational bath for which the degeneracy 
factor can be easily calculated, it still applies for cases 
with nonuniform intrabath coupling {Jj} {Jj is the in¬ 
trabath coupling strength between spin-j and spin-(j-l-l) 
in the bath). In general, the degeneracies within a certain 
magnetization sector will be lifted by introducing inho¬ 
mogeneity of {Jj}. However, the degeneracy argument 
still works since it only relies on the perturbative con¬ 
sideration at large J. As an illustration, we give in the 
Appendix an example study of a simple system having 
only = 4 bath spins, but with inhomogeneous intra¬ 
bath coupling {Ji, J 2 , da, J 4 }. 

The numerical calculations in this work are limited to 
relatively small baths with A < 16 bath spins. However, 
for the dynamics of the conventional central spin model, 
it is usually believed that A = 16 is large enough to elimi¬ 
nate finite-size effects and qualitatively similar behaviors 
with larger A can be seen. Thus, we belived that the 
above observations, namely the existence of the optimal 
J = Jopt that enhances the decoherence the most, and 
the increase of |'r(t)| when J is passing this value should 
remain for large baths with A > 16. As a result of the in¬ 
terplay between the intrabath coupling and the spin-bath 
coupling, the minimum of R{J) at Jopt occurs when J is 
comparable to {gj}. We note that the strongest dimen¬ 
sionless spin-bath coupling 2(;/(Awfluc) decreases as A 
increases, which will move Jopt/<^fluc to the left. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this work we studied the real-time quantum dynam¬ 
ics of a generalized central spin model, which extends 
the usual central spin model or Gaudin model by arrang¬ 
ing the bath spins on a circle and introducing nearest- 
neighbor homogeneous AA-type intrabath interactions. 
This model can also be viewed as an extension of the 
model studied in our previous work [36] by including the 
Ising part of the system-bath coupling. Taking advantage 
of the conservation of total magnetization and transla¬ 
tional invariance of the resulting spin bath, we could work 
in the bath’s momentum space and write down the most 
general forms of the time-dependent wavefunctions. The 
exact equations of motion for the corresponding proba¬ 
bility amplitudes in the momentum space are derived and 
integrated numerically to obtain the time evolution of the 
polarization {Sz{t)) and decoherence factor |r(t)| of the 


central spin. It can be derived from the EOMs that the 
initial decay of |r(t)| will start from at least the quadratic 
term of time t. For the case of uniform system-bath cou¬ 
pling, the EOM framework gives a transparent physical 
picture based on the momentum-conserving scattering 
processes into adjacent Lj-subsectors of the bath. 

Using the above results, we first examined the central 
spin polarization dynamics with the bath initially pre¬ 
pared in a fully polarized state. This is followed by the 
calculation of the decoherence dynamics of the central 
spin starting from an antiferromagnetic bath initial con¬ 
dition, where the high dimensions of the Hilbert spaces 
involved prevent us from dealing with large baths with 
A > 16 spins by using the EOM method. However, 
thanks to the recently developed Chebyshev expansion 
technique which proves its high efficiency in simulating 
the dynamics of the conventional central spin model [24], 
we could apply the GET approach to the current model 
to treat relatively large baths as well as ultrastrong in¬ 
trabath interactions. Gonsistence with prior works was 
found in the absence of the intrabath coupling. We re¬ 
vealed novel dynamical behaviors of the decoherence for 
finite intrabath coupling strengths. In general, the de¬ 
coherence factor \r{t)\ drops off abruptly at the initial 
stage and approaches a roughly steady value at long 
times. To illustrate this, we introduced a time-dependent 
quantity i?(t, J) to describe the time averaged value of 
|r(t)| at a given intrabath coupling strength J. Intrigu- 
ingly enough, we found that the long-time mean value 
i?( J) = i?(T'max, J) shows a non-monotonic dependence 
on J. For large enough baths with vanishing finite-size 
effect, there always exists some value of J that enhances 
the decoherence the most. However, as J increases fur¬ 
ther by passing this value, the coherence is restored and 
tends to be some fixed value that depends on the number 
of bath spins A. This dependence shows an oscillatory 
behavior as A/2 changes its parity and can be under¬ 
stood by degeneracy analysis of states in different mag¬ 
netization sectors based on the EOMs and perturbation 
theory. 

Finally, we would like to discuss possible experimental 
realizations of our model in nanoscale systems. Consider 
a modified ‘quantum corral’ [59] construction by coating 
a metallic substrate with a thin insulating layer, on the 
top of which one places iron adatoms in a ring struc¬ 
ture. The presence of the insulating layer will cut off 
the Hondo effect [60] induced by the conduction elec¬ 
trons in the metallic host, so that a direct Heisenberg 
type exchange interaction between neighboring adatoms 
might be left. By further placing an additional single iron 
adatom in the inner of the quantum corral, a model with 
Heisenberg type intrabath as well as system-bath cou¬ 
plings can be realized. The intrabath coupling strength 
J/wfluc can reach the strong coupling regime reasonably 
due to the generally larger radius of the corral compared 
to the distance between adjacent bath atoms. Although 
our model only includes an AA-type intra-ring interac¬ 
tion, we believe it is a first step towards understanding 
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this realistic system. Furthermore, the GET technique 
can, in principle, also deal with the dynamics of such a 
system with Heisenberg type intrabath coupling [61]. 

Our analysis implies that intrabath coupling could 
have significant effects on the central spin decoherence. 
The interacting central spin model and theoretical meth¬ 
ods introduced by us may stimulate further studies on the 
intrabath interaction effects in more realistic solid-state 
central spin systems. 
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APPENDIX: A 4-SPIN INHOMOGENEOUS BATH 


£;£™- 2 ) ^ ±y/2\J\, = ±|J|, and = 

E^-3) _ eigenstates with zero en¬ 
ergy E^Zi\ in the m = 2 sector, and the two 

eigenstates with zero energy E^Zp in the m = 3 sector 
mainly participate the mixing of the two subsectors. In¬ 
troducing inhomogeneitv in 1 J, I will generally lift these 



In this appendix, we will calculate the decoherence dy¬ 
namics |r(t)| for a spin bath made up of four bath spins 
with inhomogeneous intrabath couplings. The system is 
described by 

i?4 = E + 2E ■ S.- (44) 

1=1 1=1 


FIG. 5. Decoherence dynamics ]r(t)] of a central spin cou¬ 
pled to an interacting spin bath with A = 4 bath spins. 
Black: Ji = J 2 = Js = J 4 = J; green: Ji = J-\/L2, J 2 = 
J'\/0.9, J 3 = t/vDT, J 4 = blue: Ji = ./VDS, J 2 = 

J'\/0.2, Ja/tUfluc = Jv^LS, J 4 = J'\/0.2. The three sets of 
parameters are chosen such that 'Z2,. Jj is unchanged. Here 
J = Swfluc is the homogeneous intrabath coupling. 


As in Section III.B, the initial state is chosen 
as IV'(O)) = + |i))|ilil)- Diagonalizing 

the spin bath within the m = 2 (spanned by 
{[1111), jllll), jlllp, jllll), jllll), jllll)}) and m = 3 
(spanned by {jllll), jllll), jllll), jllll)}) subsectors 
gives the following eigenvalues, 

= ^^J2jf + 2a'{J,Js + J2J4), 

Eim=2) ^ _g(m=2) ^ 


degeneracies. However, if the energy splitting is of the 
order of {gj}, fast decay of jr(t)j is still possible. 

In Fig. (5), we plot jr(t)j for three sets of intrabath 
coupling configurations {Jj}, which all lie in the strong 
intrabath coupling regime in the sense that the smallest 
Jj/^fiuc ^ 2gi/a;fluc = 0.6746 for J/wfluc = 5- We keep 
Jj /4 = as a constant for each set of the param¬ 
eters. For Ji = J-\/L2, J 2 = J-\/(L9, J 3 = J\/TT, J 4 = 
JVO.S (green curve), we see that jr(f)j drops further 
compared with the homogeneous case. This can be un¬ 
derstood by looking at the eigenvalues of the spin bath: 


and 


E 


(m=3) 


2V2 


\ 



p{m=2) 
^±,1 
p(m—3) 

E±-l 


= ±7.0688 Wfluc, E. 
= ±0.3765 Wfluc, E. 


{m—2) 

±,-l 
(m=3) 
±,1 


= ±0.1789 Wfluc! and 
= ±4.9858 Wfluc- All 


where a, a' = ± 1 . 
ist two zero-energy 


We see that there 
eigenstates in the 


tor. For the homogeneous coupling Jj = J, we have 


j 2)2 _ 4 (J^ J 3 — ±2 J 4 ) 2 , the energy splittings are of the order of {gj}. In con¬ 


trast, suppression of the decay of jr(f)j is observed for 
Ji = JVts, J 2 = JV^, J 3 = JVts, J 4 = JV^ (blue 
curve). In this case, the spectrum in the m = 2 sector 
does not change as Ji J 3 ± J 2 J 4 = 2J^. However, Ep}Zp 
raise up to ±2.2361 Wfluc ^ gj, which consequently sup¬ 
press the mixing of the two subsectors. 


(46) 


always ex- 
m = 2 sec- 









